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Abstract 

We revise the problem first addressed by Braunstein and co-workers 
(Phys. Rev. Lett. 83 (5) (1999) 1054) concerning the separability of 
very noisy mixed states represented by general density matrices with 
the form p t = (1 — e)Md + ep\. From a detailed numerical analysis, it is 
shown that: (1) there exist infinite values in the interval taken for the 
density matrix expansion coefficients, —1 < c aii ... jajv < 1, which give 
rise to non-physical density matrices, with trace equal to 1, but at least 
one negative eigenvalue; (2) there exist entangled matrices outside the 
predicted entanglement region, and (3) there exist separable matrices 
inside the same region. It is also shown that the lower and upper 
bounds of e depend on the coefficients of the expansion of p\ in the 
Pauli basis. If p\ is hermitian with trace equal to 1, but is allowed to 
have negative eigenvalues, it is shown that p e can be entangled, even 
for two qubits. 



The literature of Nuclear Magnetic Resonance quantum computing has 
been one of the most fruitful for the past few years, since the discovery of 
pseudo-pure states by Gershenfeld and Chuang [1] and Cory et al. [2] . Exper- 
iments performed by different groups reported various implementations of al- 
gorithms [3, 4, 5], simulations [6] and quantum entanglement [7, 9], including 
teleportation [8] and other quantum effects [10, 11]. These results contrast 
with the landmark paper published by Braunstein and co-workers in 1999 
[12] in which it is argued that, with the present stage of NMR technology, 
all those experiments could be interpreted through usual classical correla- 
tions between spins. These conclusions were further extended by Linden and 
Pospescu [13]. The latest of such experiments (to the best of the authors 
knowledge) was performed in 2003 by Mehring et. al [14] who reported en- 
tanglement of nuclear and electron spins in a molecular single-crystal, using 
the combined techniques of NMR and EPR. 

Reference [12] considered arbitrary density matrices for iV qubits in the 
form 

p e = (1 - e)M d + e Pl (1) 

where d = 2 N is the dimension of the Hilbert space for N qubits and Md = 
I/d the maximally mixed density matrix. I is the identity matrix in the 
space of N qubits, and p\ an arbitrary density matrix. Matrices of this form 
are expanded in a basis of Pauli density matrices: 

P = ^c ai ... aN a ai <g> ... <g> a aN (2) 
where {a s } = {0,i s } = {0,1,2,3} and "s" indicates the s— th qubit with 
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the sum made over repeated indices. Normalization imposes C0...0 = 1 and 
the other coefficients are in the interval — 1 < c ai _ aN < 1. From this, after 
a transformation for an overcomplete basis, it is established that, for the 
case N = 2, taking the minimum value of the coefficients, c aijQ2 = — 1, the 
bound e < 1/15 limits the region below which p e is separable. Generalization 
for arbitrary N leads to e < 1/A N . Since typically e w 10~ 5 in NMR room 
temperature liquid state experiments, the conclusion is that so far no entan- 
glement has ever taken place in NMR experiments, a conclusion which has 
been revised by others [15, 16]. One important observation is that in Ref. 
[12] this bound is assumed to hold independently of p\. 

However, since the density matrix in Eq. (1) is arbitrary we can ap- 
ply the reasoning to the simplest case: N = 1 (of course, this involves no 
entanglement!). Taking c ai = —1, this leads to: 

pN =^2\-l-i 2 ) ® 
We see that this matrix satisfies the condition Tr(p) = 1, but its eigenvalues 
are Ai = anc ^ ^ 2 = Therefore, A 2 < 0, and Eq. (3) cannot rep- 

resent a density matrix of a physical system [17]. The same can be observed 
for the cases N = 2, with c QliQ2 = — 1, and N = 3, with c ait0l2t0l3 = — 1, whose 
respective density matrices are: 



I -2 -2 + 2i -2 + 2i 2i \ 



1 

Pn=2 — 7 



(4) 



—2 — 2i 2 -2 

-1-1% -2 2 

V -2% 2 / 

also satisfying the normalization constraint Tr(p) = 1, but with eigenvalues 

Ai = =^,A 2 = =2^ < 0,A 3 = A 4 = 1, and 
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X 



1 

pN=3 = gX 

/ -6 -4 + Ai -4 + Ai Ai -4 + 4i Ai Ai 2 + 2i \ 

—A — Ai 2 -4 -4 -2 + 2i 

-4-4i -4 2 -4 -2 + 2i 

-Ai 2 —2 — 2i 

-A-Ai -4 -4 -2 + 2i 2 

-Ai -2-2i 2 

-Ai -2-2i 2 

V 2 - 2i 2 

(5) 

with Tr(p) = 1, and eigenvalues Ai = ~ 8 + 6 ^ , A 2 = ~ 8 ~ 6v ^ < 0, A 3 = A 4 = 
A5 = 4+ g^ , A6 = A7 = As = 4 ~g^ . Since a physical density matrix must 
be a positive operator [18], none of the above matrices can be regarded as 
representing a physical system. 

In order to find valid intervals, consider an example where all the coeffi- 
cients are equal to some constant c ai) ... )Q , JV = c in Eq.(2), and let us impose 
A > for the eigenvalues of the resulting matrix. This leads to the following 
intervals: for N = 1, -0.58 < c < 0.58; for JV = 2, -0.15 < c < 0.33; for 
N — 3, —0.05 < c < 0.15, and so on. This case defines only one possible set 
of values for the coefficients, but will be useful to derive a number of results. 
Note that, for = 2 and c QljQ2 = c, the intervals —1 < c < —0.15 and 
0.33 < c < 1 define infinite non-physical matrices within — 1 < c < 1. 

With the new intervals, the eigenvalues for the case N = 2, for c = —0.15, 
are all positive and p satisfies tr(p) = 1: 

{Xi} = {0.007596,0.267404,0.362499,0.362500} 
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The same is true for N = 3 and c = —0.05: 

{Ai} = {0.003798,0.133702,0.165401,0.165401,0.165401, 

0.122099, 0.122099, 0.122099} 

Note that, except for N = 1, the intervals are asymetric. If we write 
the intervals in the form — I/An < c Ql ... ajv < 1/Bn (where, for instance, 
{A 2 ,B 2 ) = (6.67,3.03), {A 3 ,B 3 ) = (20.00,6.67)), applying the procedure 
of [12] we find that the lower bound for e in the case N = 2 would be 
e < A 2 /15 = 0.44, larger than the previous e < 1/15 and, in the case N = 3, 
e < ^3/63 = 0.32, also larger than the previous e < 1/64. According to Ref. 
[12], for N = 2, a matrix p e with e > 0.33 would be entangled, but according 
to the above there exist matrices in this interval which are separable. It is 
important to recall that in that reference, the bounds are independent of p±. 
This is not a problem of basis choice, for it also appears in the second case 
considered in Ref. [12] , in which a continuous overcomplete basis is used. In 
this case, from Eq. (9) of [12] one can considers those w(fii, . . . , Un) which, 
for a > 1, satisfy the relation u>(ni, . . . ,njv) > — 2 2N ~ 1 /a(4n) N , we see that 
p e is sepparable for e < a/(a + 2 2N ~ 1 ). In particular, for N = 2 and a = 6 we 
see that p t will be separable for e < 3/7, but according to Ref. [12], matrices 
p t with e > 0.33 are entangled. 

Next, let us present an explicit example of separable p t forN = 2 within 
the entangled region given in Ref. [12]. Let pi be a matrix which, in the 
Pauli basis, is defined by the coefficients c aiiQ , 2 = c = —0.15. Thus, using Eq. 
(2), one has: 



5 



/ 0.1375 -0.0750 + 0.0750i -0.0750 + 0.0750i 0.0750i \ 



-0.0750 - 0.0750i 0.2875 -0.0750 

Pl = -0.0750 - 0.0750i -0.0750 0.2875 

V -0.0750i 0.2875 / 

(6) 

This matrix satisfies Tr(p) = 1 and has eigenvalues Ai = 0.0076, A 2 = 



0.2674, A 3 = A 4 = 0.3625. Therefore, it is a true density matrix. Now, if 
we take e = 0.40 (that is, entangled according to Ref. [12]) and replace in 
Eq. (1), one obtains 

/ 0.2050 -0.0300 + 0.0300i -0.0300 + 0.0300* 0.0300i \ 



-0.0300 - 0.0300i 0.2650 -0.0300 

-0.0300 - 0.0300i -0.0300 0.2650 

V -0.0300i 0.2650 / 

(7) 



which is also a true density matrix, since Tr(p e ) = 1, and has positive eigen- 
values: Ai = 0.1530, A 2 = 0.2569, A 3 = A 4 = 0.2950. But, according to Peres 
criterium [19], this matrix is separable. To show that, let us expand p e in 
the following basis of density matrices: 

pi = l{ I i ) ; p2 = l{) 7 ) ; p * = I ( 1 1 ) ; p4 = ( I o ) 

(8) 

After partial transposing [19], one obtains the following eigenvalues: Ai = 
0.1530, A 2 = 0.2569, A 3 = A 4 = 0.2950. Since for two qubits the Peres cri- 
terium is a necessary and sufficient condition for separability [20], we con- 
clude that p e is separable. The same result can be reached using the same 
procedure given in Ref. [12]. To see that, we need the coefficients of p t in the 
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Pauli basis, d aij0 , 2 , which are equal to ce for the case the coeffcients of p\ are 
all equal to c (see the Appendix). Replacing into Eq. (4) of Ref. [12], one 
has: 

1 0.15e 0.15e 
WiWi + d i0 Wj + Wid 0j + dij = 0.15e > 

from which one obtains the separability condition: e < 0.44, satisfied by the 
previous p t , with e = 0.40. This is not the only case; there exist infinite 
separable matrices p t with e > 0.33. This occurs because the upper and 
lower bounds are not fixed, as in Ref. [12], but change with p\. This is why 
in Eq.(9) e > 0.44 must be considered, and not e > 0.33, as the correct upper 
bound. 

Let us now build a matrix p 1 with the following set of unequal coefficients: 
coi = -0.07; c 02 = -0.07; c 03 = -0.07; c w = -0.06; c n = -0.83; c 12 = 
-0.03; ci 3 = -0.03; c 20 = -0.06; c 21 = -0.03; c 22 = -0.03; c 23 = 
—0.03; c 30 = —0.06; c 3 i = —0.03; c 32 = —0.03; c 33 = —0.03 which has eigen- 
values {\i = 0.000996,0.078152, 0.448038,0.472814}. Thus, for e = 0.13 one 
obtains: 



Pe = 



( +0.2448 -0.0032 + 0.0032i -0.0029 + 0.0029i 

-0.0032 - 0.0032i +0.2513 -0.0279 

-0.0029 - 0.0029i -0.0279 +0.2506 

\ -0.0260 - 0.0019i -0.0010 - O.OOlOi -0.0013 - 0.0013i 



(9) 



-0.0260 + 0.0019i \ 
-0.0010 + 0.0010i 
-0.0013 + 0.0013i 
+0.2532 J 

with eigenvalues: Ai = 0.217629, A 2 = 0.227660, A 3 = 0.275745, A 4 = 
0.278966. Applying Eq. (4) of [12], one has: 

w\.wi + dioWi + widoi + dn = —0.0218 < 



where dij are the coefficients of p € in the Pauli basis (see Appendix). There- 
fore, p e is entangled, but with e outside the bounds established in Ref. [12]. 
The appropriate procedure to determine whether NMR can or cannot im- 
plement entanglement consist in finding the matrix pi, for instance from 
quantum state tomography experiments [18], determine the expansion co- 
efficients of Eq.(2) and, from those, derive the respective lower and upper 
bounds which will be valid only for a particular set of matrices p e . 

Finally, we will show that with the only two requirements of p\ being 
hermitian and satisfying Tr(p 1 ) = 1, a matrix p t of two qubits can be en- 
tangled, for usual values of e. In order p t to be a true density matrix the 
only requirement on p\ is that it must be hermitian and have trace equal 1. 
Imposing non-negative eigenvalues (EV) for p e , one has: 

EV{ Pe } = EV{(l-e)M d + e Pl }>0 (10) 

Since [M d + p 1} p x ] = 0, [M d + p u M d ] =0 and [pi, M d ] = 0, the EV of the 
sum is equal to the sum of EV: 

EV{ Pl } > -i^ (11) 

That is, for the usual values e ps 1CT 5 , the absolute values of the EV of pi 
can actually be very large. 

Let us restrict ourselves to the case N = 2 and pick the following partic- 
ular set of values for the coefficients: c Ql . jQ2 = —666.66 (note that since here 
P i does not have to be a density matrix, these are valid coefficients). From 
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this, we obtain: 

/ -499.745 -333.33 + 333.33i -333.33 + 333.33i 333.33i \ 



-333.33 - 333.33i 166.915 -333.33 

-333.33 - 333.33i -333.33 166.915 

V -333.33i 166.915 / 

(12) 



which has eigenvalues Ai = -1077.089496, A 2 = 77.599495, A 3 = 500.244999, 
A 4 = 500.245000. From this, and setting e = 0.0002 in Eq.(l) we find the 
matrix: 

/ 0.15 -0.0666 + 0.0666i -0.0666 + 0.0666i 0.0667i \ 



-0.0666 - 0.0666i 0.2833 -0.0666 

-0.0666 - 0.0666i -0.0666 0.2833 

V -0.0667i 0.2833 / 

(13) 



which has eiganvalues A x = 0.034532, A 2 = 0.265469, A 3 =0.3499989,A 4 = 
0.3499990, and therefore is a good density matrix. According to the criterium 
of Ref. [12] this matrix would be entangled. 

Summarizing, the method of Ref. [12] to decide whether NMR can or 
cannot produce entanglement has been revised. A numerical analysis has 
shown that the interval ±1 for the expansion coefficients is not generally 
correct, and that the bounds for entanglement are dependent on the matrix 
p\. The non observation of this fact leads to the following problems: (1) 
existence of density matrices with negative eigenvalues; (2) existence of en- 
tangled matrices outside the predicted entanglement region; (3) existence of 
separable matrices inside the same entanglement region. We also show that 
with the requirements of p\ being hermitian with trace equal 1, but allowed 
to have negative eigenvalues (a possible experimental situation, as reported 
in Ref. [21]), p t can be entangled, even for two-qubits. However, we were not 
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able to produce entangled two qubits p e , for p\ being a true density matrix 
(that is, a positive operator), what is in accordance with one conclusion of 
Ref. [12]. Cases for higher numbers of qubits must be analyzed in separate, 
and they may lead to further insights into this problem. 
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1 Appendix 

Any given two-qubit matrix p possesses the following elements in the expan- 
sion in the Pauli basis: 

Pn = i(l + c 03 + c 3 o + C33) P12 = i(c 01 - ic 02 + c 3 i - ic 32 ) 

Pl3 = i( Cl + C 13 - ic 20 - ^C 23 ) Pl4 = \{cu - tC 12 - ic 21 ~ C 22 ) 

P21 = i( C 01 + ic 02 + c 31 + IC32) P22 = \(l - C 3 + C 30 - C33) 

P23 = |(cn + ic 12 - ic 21 + C 22 ) P24 = |(C1 - C13 - ic 20 + «c 23 ) 

P31 = i( C 10 + C 13 + «C 20 + «C 23 ) P32 = \(CU ~ IC 12 + ic 21 + C 22 ) 

P33 = j(l + c 03 - c 30 - c 33 ) p 34 = i(c i - ic 02 - c 3 i + ic 32 ) 

P41 = i(cn + iCi2 + ic 2 i - C 22 ) p 42 = |(ci - C13 + ic 20 - ic 23 ) 
P43 = i( C 01 + ^ C 02 - c 31 - ^C 32 ) p 44 = j(l - C 03 - C 30 + C33) 

Given pi for iV = 2 and given e, the matrix p t will be defined by Eq. 
(1), which can in turn also be expanded in the Pauli basis. Let d au0lj be the 
expansion coefficients of p 6 . Thus, for instance, for the element (1, 1) we find 
the following relation: 

+ |(1 + C03 + C 30 + C33) = ^(1 + ^03 + ^30 + d 33 ) (14) 

Letting d aija2 = ec ai ^ a2 the above equation will be valid for any value of e, 
particularly for c 03 = c 30 = c 33 = c and do3 = d 3 o = d 33 = d. It is easy to 
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show that by setting c au a 2 = c e d ai:0 , 2 = d = ce one obtains the same result 
for any two corresponding elements. 
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